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REFLECTION  OF  SPHERICAL  WAVES  AT  WEAK  INTERFACES 


Zhurnal  Tekhnicheskoy  Fizlki  L.  Brekhovskikh 

(Journal  of  Technical  Physics) 

Vol  XVIII,  No.  4,  1968, 

Pages  473-483 


The  field  of  an  acousclc  or  radiowave  point  source 
at  a  flat  interface  is  investigated  in  the  case  where  the 
properties  of  the  media  on  both  sides  of  the  interface  are 
similar.  The  limits  of  applicability  of  the  existing  theories 
are  pointed  out.  Criteria  are  found  for  which  transition 
layers  may  be  replaced  by  interfaces. 

Introduction 

We  encounter  the  problem  of  reflection  of  spherical  waves  when 
investigating  the  field  of  an  electromagnetic  or  acoustic  point  source 
in  the  presence  of  an  interface  between  two  media.  This  problem  is  the 
subject  of  a  large  number  of  works  in  which  It  has  been  investigated 
basically  from  two  points  of  view:  a)  investigation  of  the  field  at 
strongly  reflecting  interfaces, which  in  the  case  of  electromagnetic  waves 
corresponds  to  propagation  of  radiowaves  over  good  conducting  surfaces 
[1-4] ;  b)  investigation  of  the  field  in  the  wave  range  with  any  proper¬ 
ties  of  the  media  (5,  6,  13].  These  two  cases  actually  overlap  each 
other.  The  cited  papers  as  a  group  represent  a  complete  Investigation 
of  the  problem  satisfying  the  majority  of  practical  and  experimental 
questions. 

However,  returning  to  the  solutions  obtained  in  these  papers,  let 
us  note  that  they  give  infinitely  increasing  values  for  the  amplitudes 
of  the  reflected  waves  on  approximation  of  the  index  of  refraction  n  to 
one  at  the  interface  (see  part  1  in  [2]  and  also  [5,  6,  13]).  Clarifica¬ 
tion  of  the  nature  of  the  anomaly  obtained  here,  determination  of  the 


i 


limit*  of  applicability  of  tha  existing  theories  and  ;derivation  of  formu¬ 
las  which  are  valid  for  n  aa  cloae  as  deaired  to  'me  —  all  are  of  theo¬ 
retical  end  practical  interest*  Zn  particular,  when  investigating  radio 
wave  propagation  in  the  atmosphere  frequently  it  is  necessary  to  deal 
with  their  reflection  from  layers  which  may  be  reduced  to  interfaces  be¬ 
tween  media  with  similar  dielectric  constants.  The  index  of  refraction 
for  these  interfaces  may  differ  from  one  by  only  small  fractions  of  a 
%  [7]*  We  encounter  the  taste,  case  on  reflection  of  sound  from  layers  of 
discontinuity  in  the  sea.  As  is  known,  the  latter  are  comparatively 
thin  layers  extending  in  the  horizontal  direction  with  large  temperature 
and  salinity  gradients.  The  index  of  refraction  equal  to  the  ratio  of 
the  propagation  rate  of  sound  in  the  water  above  and  Sbelow  a  layer  of 
discontinuity  also  differs  from  one  by  only  a  few  percent  or  fraction  of 
a  percent. 

The  interfaces  of  media  with  similar  properties  and,  consequently, 
an  index  of  refraction  close  to  one  will  be  called  weak  interfaces  in 
the  remainder  of  this  article. 

Let  us  note  that  the  formulas  obtained  by  Ryazin  [12]  for  the  Hertzian 
vector  are  also  valid  for  weak  interfaces.  However,  they  refer  only  to 
the  case  >vhere  the  transmitter  and  receiver  are  both  at  the  interface. 

§  1.  Limits  of  Applicability  of  the  Existing  Theories 

Apparently  we  were  the  first  to  note  [13]  the  anomaly  obtained  when 
n  1,  and  also  the  fact  that  investigation  of  the  field  of  a  reflected 
wave  by  expansion  in  a  series  with  respect  to  powers  of  l/fc^tf  (which  is 

done  by  the  saddle  point  method,  used  in  references  [1,  5,  6])  is  possible 
for  all  glancing  angles  only  on  fulfillment  of  the  condition  ([13],  equa¬ 
tion  (28)  : 


where  R  is  the  distance  from  the  image  radiator  to  the  receiver.  If  the 
latter  condition  is  not  fulfilled,  the  saddle  point  method  gives'  correct 
values  for  the  fields  only  at  sufficiently  large  glancing  angles  x 
(figure  1)  satisfying  the  condition 

1.  (2) 

The  inapplicability  of  the  saddle  point  method  for  n  cl&je  to  one 
and  small  glancing  angles  may  be  understood  frost  some  very  obvious  argu¬ 
ments  which  we  consider  it  useful  to  present.  On  reflection  of  e  spheri¬ 
cal  wave. all  points  of  the  interface  do  not  play  an  essential  role,  but 
only  some  effective  zone  [14]  in  the  shape  of  an  ellipse  the  area  of  which 


approaches  zero  on  conversion  to  geometric  optics  where  it  is  possible  to 
talk  about  reflection  of  the  wave  from  a  certain  point  on  an  interface. 

In  accord  v/ith  this  the  field  of  the  reflected  wave  at  an  arbitrary  point 
in  space  is  the  result  of  waves  going  only  in  the  directions  which  corre¬ 
spond  to  straight  lines  connecting  this  point  with  all  parts  of  the  effec¬ 
tive  zone.  These  will  be  the  directions  deviating  by  small  angles  on  the 
order  of  I/Sk^R  and  less  from  the  direction  of  the  beam  reflected  by  the 

laws  of  geometric  optics.  Application  of  the  saddle  point  method  is  pos¬ 
sible  if  the  reflection  coefficient  and  its  derivative  with  respect  to  the 
angle  may  be  considered  to  have  low  variation  in  this  range  of  angles . 

However,  from  the  known  expression  for  the  reflection  coefficient 
of  plane  waves  reflected  from  the  interface  of  two  media 

B  (X)=  (3) 

m  sin  y  \V-  —  cos-  y, 

we  find  that  for  n  close  to  one  and  small  x  the  first  and  second  deriva¬ 
tives  of  B(x)  way  be  as  large  as  desired.  Here  m  *  in  electrodynamics 

and  m  *  Pj/Pq  (£lie  density  ratio  of  the  media)  in  acoustics. 

Application  of  the  saddle  point  method  gives  geometric  optics 
(acoustics)  in  the  first  approximation  and  the  corrections  to  it  in 
subsequent  expressions.  For  n  close  to  one,  fulfillment  of  condition 
(2)  ensures  that  the  corrections  will  be  small.  Therefore,  it  is  also 
a  condition  of  applicability  of  geometric  optics  in  the  description  of 
the  field  of  a  reflected  wave. 

§  2.  The  Field  of  a  Point  Source  Located  at  a  Weak  Interface 

Let  us  investigate  the  field  of  a  point  source  located  at  a 
height  of  Zq  (which  may  also  be  zero)  of  a  flat  interface  at  the  point 

(v,  2 ) .  In  the  electromagnetic  case  a  vertical  dipole  will  be  used  as 
the  radiator,  and  in  the  acoustic  case,  a  pulsating  sphere  of  infinitely 
small  radius.  The  field  will  be  characterized  by  a  scalar  function  tf 
which  is  the  vertical  component  of  the  Hertzian  vector  or  the  acoustic 
potential  respectively. 

It  may  be  demonstrated  (9]  that  for  any  relationship  of  fixed  media 
we  will  have 


/f| 


(4) 


for  the  field  in  the  upper  medium  where  the  first  term  is  the  direct  wave 
and  the  second  is  the  reflected  wave.  Here 


f  mio  —  lt 
J  m6o-+~6j 


'  (5) 


Figure  1.  R ^  —  the  distance  from  the  reception  point  (r,z) 

to  the  radiator  Q;  R  —  the  distance  from  the  same  point  to 
the  image  radiator  x  *s  the  glancing  angle  formed  by  the 
path  of  the  beam  with  the  Interface. 

Here  '60==V£i  — V;  bl=^~kli;  kQ  and  kj  are  the  wave  numbers  res¬ 
pectively  in  the  upper  and  lower  media,  and  m  has  the  same  value  as  in 
(3),  and  C  -  z  +  aQ. 

Our  problem  is  to  investigate  expression  (5)  for  n  »  k^/k^  close 

to  one. 


Let  us  expand  the  expression  under  the  integral  sign  in  (5)  in  a 
series  with  respect  to  powers  of  n2  -  1.  For  this  purpose  we  shall 
introduce  the  notation 


Then 


V  V  * 

(6) 

•  *  CO 

mbo  —  bt _ m  —  V  n  , 

«= 0 

(7) 

Substituting  this  expression  in  (S)  and  considering  (6),  for  the 
reflected  wave1  we  obtain 

^The  radiuB  of  convergence  of  the  series  (7)  is  equal  to  one.  There¬ 
fore,  representation  of  <p  in  the  form  of  (8)  is  possible  only  in  the  case 
where  the  integration  patK  is  selected  in  (5)  on  which  |a;|  <1  everywhere, 
that  is,  \b0\z  >  | k^1  -  k-2 I •  For  example,  we  can  take  for  this  path  the 
path  which  runs  from  £  *  6  first  along  the  real  axis,  then  by  way  of  the 
point  £  *  kg  along  the  semicircle  of  radius  greater  than  |/52T2-'  757*]* 
lying  in  the  fourth  quadrant;  and  then  back  to  the  real  axis./7  Since  the 
expression  under  the  integral  sign  in  (5)  has  nc  singularities  there  are 
no  obstacles  to  converting  to  this  path  of  integration.  It  is  assumed  that 
in  the  integrals  (9)  obtained  after  expansion  in  the  sferies,  the  integration 
path  again  coincides  with  the  real  axis  from  which,  of  course,  the  values 
of  the  integrals  are  not  changed. 


(8> 
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(16) 


)  dt—  j  —j^—tdt. 


Here  the  second  integral  is  taken  directly,  and  on  substitution  of 
the  limits  it  gives  /  (l  r) 

%{e  e  (17) 

To  calculate  the  first  integral  let  us  replace  the  variables  by 
the  formula 

*o(^i  —  r)=at*,,  from  which  ■^^•=2xrfa,  • 


*.<A-r) 


f  i__cV=2e'*‘'  f 


We  shall  assume  that 


S«i. 


and  then  on  the  entire  integration  path  ors~<^>5:0r,  as  a  result  of 
which  the  first  term  under  the  square  root  in  (18)  may  be  neglected. 


As  a  result  we  obtain 


where 


and  C(u )  and  S(u )  are  the  Fresnel  integrals.  Substituting  this  result 
in  (14)  and  considering  (19)  we  obtain  the  expression  for  I-.  On  sub¬ 
stituting  this  in  the  expression  for  the  reflected  wave  (13;  we  have 

tif!'  (m  —  1  im(1  —  rz-)kltr  ,  —,T I r>  .  -p  l-**/)"]]  /on 

?,  =  —  |  - [1  -«««  (C+tS-  -2-  JJp  (->:) 

2 

In  electrodynamics  where  m  -  n  ,  the  first  term  in  the  braces  is 
k^r  times  less  than  the  second  and  consequently  may  be  omitted.  In 
acoustics  we  are  interested  only  in  the  cases  where  m  -  1  is  on  the  order 
of  or  less  than  n2  -  1  so  that  the  first  term  may  again  be  omitted.* 

As  a  result  the  reflected  wave  is  written  as 


n  +  1. 


len  m  -  1  is  not  awall,  we  do  not  obtain  a  singularity  when 
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where 


9r 


B 


e 


ttj  /: 


R 


(22) 


n _ »n  —  1  Ri  r 


1  — stius 


(&) 


and  in  electrodynamics  the  factor  mj(m  +  l)J«n!/(«1  + 1)2  may  be  replaced 
by  1/4. 


§  3.  Discussion  of  the  Results 


Expression  (22)  represents  an  ordinary  spherical  wave  reflected  with 
a  reflection  coefficient  of  B.  It  may  be  shown  that  for  u2  »  1  it  co¬ 
incides  with  the  reflection  coefficient  of  plane  waves  (3)  if  the  latter 
is  also  expanded  inaseries  with  respect  to  powers  of  nz  -  1  and  is 
limited  to  the  first  power.  This  must  be  expected  inasmuch  as  the  con¬ 
dition  u2  »  1  which  may  also  be  written  in  the  form 

(24) 

is  equivalent  to  condition  (2)  on  fulfillment  of  which  geometric  optics 
is  valid. 

Let  us  introduce  the  new  variable 


x  =  T-ut 


in  place  of  u,  and  let  us  use  the  asymptotic  expansion  of  the  Fresnel 
integrals  [8J  with  respect  to  powers  of  l/«. 

Then  we  obtain 


C-*-iS- 


1  H-i 


V2 


1  l ’• 
2=7  (  a 


1-3 


2x  ~(2xf  1 ' 


"•)* 


Substituting  this  in  (23)  and  considering  that  sin  x  **  C/ff,  drop¬ 
ping  small  values  of 


d  _  m(l-n-) 

— (m 


(25) 


we  obtain  for  B  the  same  expression  as  ve  obtained  on  expanding  (3)  with 
respect  to  powers  of  nz  -  I  limiting  ourselves  to  the  first  term  and 

dropping  the  term  2LH.I  in  accord  with  the  above  indicated  arguments. 

tn  -i- 1 

When  u  «  1  in  (23)  all  terms  in  brackets,  except  the  1,  may  be 
neglected  as  a  result  of  which  the  reflection  coefficient  B  again  assumes 
a  simple  form. 
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It  is  of  interest  to  trace  the  value  of  the  ratio 
values  of  w. 


for  all 


According  to  (23)  and  (25) ,  we  have 

^  =  —  wu2  j^l — viua  1  (cL\~iS— (26) 

The  modulus  of  this  quantity  is  depicted  in  figure  2  as  a  function 
of  u  where  the  logarithmic  scale  is  taken  along  the  axes.  With  an  increase 
in  u,  j S/B  I  increases  continuously.  This  increase  first  proceeds 
according  to  the  law  u *•  which  on  the  assumed  scale  corresponds  to  a 
straight  line.  For  large  u,  |b/Bq|  approaches  zero  which  conditions  the 
applicability  of  geometric  optics,  that  is,  ooinddano*  of  the  reflection 
coefficients  of  spherical  and  plan*  waves.  Thi3  result  refers  to  the 
case  of  a  sufficiently  thin  transition  layer  (see  §  4)  when  the  latter 
may  be  replaced  by  an  interface.  For  the  case  of  thick  transition  layers 
we  have  only  the  calculations  of  reflection  of  plane  waves. 

We  may  consider  that  the  limits  of  applicability  of  these  calcula¬ 
tions  to  spherical  waves  are  also  determined  by  the  condition  (24). 


Figure  2.  Modulus  of  the  ratio  #/£«  as  a  function 
of  a  parameter  u.  B  is  the  reflection  coefficient 
•nsui&g  from  the  exact  theory  of  reflection  of  a 
spherical  wave;  Bq  is  the  reflection  coefficient 
for  a  plane  wave  (geometric  optics  or  acoustics). 

Let  us  present  some  examples  where  condition  (24)  is  not  fulfilled. 
For  the  case  of  radio  waves  with  location  of  the  weak  interface  at  a 
height  of  400  meters  above  the  earth  and  at  distances  of  50  km  the  value 
of  k^~tR  is  10  when  X  «  2  meters  and  2  when  \  ■  10  meters. 


For  a  detailed  bibliography  see  the  article  by  Vvedenskiy  and 
Arenberg  [ 7  j . 
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In  the  latter  case  the  field  of  the  reflected  wave  must  be  calcu¬ 
lated  by  formulas  (22)  and  (23)  obtained  above. 

In  hydroacoustics  with  a  depth  of  the  layer  of  discontinuity  of 
30  meters  and  at  a  distance  ofi?  =  2  km  the  reflected  wave  cannot  be 
calculated  by  the  laws  of  geometric  acoustics  f or  X  >  1  meter. 

In  practical  respects  investigation  of  the  propagation  of  sound 
or  radio  waves  in  a  layer  bounded  on  one  side  by  a  weak  interface  is 
of  interest.  We  demonstrated  earlier  [9]  that  in  the  most  general  case 
the  field  in  the  layer  may  be  broken  down  into  an  infinite  sum  of  waves 
reaching  the  receiver  after  a  different  number  of  reflections  from  the 
interface.  In  the  investigated  case  the  amplitude  of  a  wave  undergoing 
k  reflections  from  a  weak  interface  will  have  a  factor  of  in  ~  1)* 
from  whifch  it  follows  that  in  our  approximation  it  is  necessary  to  consider 
only  the  following  six  waves  not  having  one  or  only  one  reflection  from 
this  interfaces  1)  a  direct  wave,  2)  a  wave  reflected  from  an  ordinary 
interface,  c)  a  wave  reflected  from  a  weak  interface.  A)  a  wave  reflected 
from  an  ordinary  and  then  a  weak  interface,  5)  the  same  with  the  inverse 
order  of  reflections,  6)  a  wave  reflected  twice  from  an  ordinary  inter¬ 
face  and  once  from  a  weak  interface.  The  second  of  the  enumerated  waves 
is  calculated  by  the  formulas  obtained  above  [13],  and  we  investigated 
the  third  one  above.  The  fourth,  fifth  and  sixth  reduce,  to  the  third 
if  the  ordinary  interface  is  an  absolutely  reflecting  interface  as,  for 
example,  the  surface  of  the  water  in  the  hydroacoustic  case.  Here  the 
wave  amplitude  will  be  calculated  by  formulas  (22)  and  (23)  where  formula 
(20)  will  be  valid  for  u  if  by  £  we  mean  the  projection  qf  the  path  followed 
by  the  corresponding  beam  to  the  3-axis,  and  R  «  +  £  .  In  addition, 

if  for  an  absolutely  reflecting  interface  the  reflection  coefficient  Bq  m 
-1,  as  in  the  above  indicated  hydroacoustic  case,  when  the  second  inter¬ 
face  is  not  absolutely  reflecting,  the  problem  is  somewhat  more  compli¬ 
cated,  but  it  may  also  be  solved  with  the  help  of  expansions  of  the 
solution  in  a  series  with  respect  to  powers  of  n  -  1  and  subsequent 
application  of  the  saddle  point  method.  However,  we  shall  not  discuss 
this. 


Replacement  of  the  Transition  Layer  by  an  Interface 


In  the  majority  of  practical  cases  a  transition  layer  occurs  in 
place  of  a  weak  interface.  The  question  arises  as  to  the  cases  in  which 
this  layer  may  be  idealized  as  an  interface.  It  is  known  that  in  ordin¬ 
ary  cases  this  may  be  found  if  the  thickness  of  the  layer  is  small  by 
comparison  with  the  wavelength.  Weak  interfaces  also  have  a  peculiarity 
in  this  ca3e  since  the  corresponding  condition  for  them  turns  out  to  be 
weaker.  To  derive  the  latter  let  us  use  some  obvious  arguaents. 


Let  the  transition  layer  with  a  thickness  l  be  parallel  to  the 
plane  3=0.  Let  us  consider  reflection  of  the  plane  wave  incident  on  it 
at  a  glancing  angle  x*  We  are  interested  in  the  caaes  where  l  is  com¬ 
parable  to  the  wave  length  X  or  greater  than  it  where  the  difference  in 
propagation  rates  on  the  upper  and  lower  boundaries  of  the  layer  is 
small.  Under  these  conditions  the  propagation  of  a  plane  wave  will  be 
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subject  to  the  laws  of  geometric  optics  so  that  cue  dependency  of  the 

phase  of  the  wave  on  the  coordinate  z  will  be  given  by  the  factor  exp 

[i  f  k  da],  The  extent  of  the  layer  in  the  a  direction  may  be  neglected 
z 

if  the  phase  inroad  in  the  thickness  of  the  layer  is  snail,  That  is,  if 

i 

f  k,dz<  1. 

0 

The  latter  condition  will  be  reinforced  by  replacing  k by  its 
maximum  value  in  the  layer  as  a  result  of  which  it  is  written  as 

(fc  )mnv  l  «  1.  (27) 

z  max 

With  a  monotonic  change  in  k&  within  the  boundaries  of  the  layer 

(only  this  case  is  of  interest  to  us)  the  maximum  value  of  k  is  reached 

z 

on  one  boundary  of  the  layer  as  a  result  of  which  the  preceding  condition . 
is  entirely  equivalent  to  the  following  two  conditions 

(Vo 1  K<  1  CW  <<  1»  (28) 

where  the  indices  0  and  1  refer  to  the  media  located  above  and  below  the 
layer  respectively.  Since  (&  )^  *  27r/XQ  sin  x  and  ( k :  )^  ■  2v/\^  sin  x^ 

where  x^  the  glancing  angle  of  the  refracted  wave,  condition  (28)  is 
rewritten  in  the  form 

/ sin  y.  >0/277,  /sin  yA  X}j2~.  (29) 

In  the  esse  of  weak  interfaces  X^  *  X*  Xq  "  and  then  both  condi¬ 


tions  reduce  to  one  condition 


/  sin/.<^X/2". 


For  small  glancing  angles  the  last  condition  may  be  fulfilled  for  t 
comparable  to  and  even  greater  than  X. 

The  conversion  from  plane  waves  to  spherical  waves  causes  complica¬ 
tions  at  first  glance.  Actually,  on  expansion  of  the  spherical  wave  in 
plane  waves  ((2] ,  and  also  [3],  page  943)  the  wave  must  be  considered  with 
a  wave  vector  component  along  the  z  axis  as  large  as  desired  and  Imaginary. 

For  large  k  the  condition  (27)  is  not  fulfilled.  However,  this  difficulty 
z 


-  10  - 


disappears  if  we  consider  the  field  in  the  wave  zone  (fc ^  »  1).  Then 

as  was  pointed  out  in  §  2,  only  plane  waves  will  play  a  significant  role 
which  have  glancing  angles  close  to  the  angle  x  depicted  in  figure  1.  This 
angle  must  be  substituted  in  condition  (30) . 

Appendix 

The  integrals  (9)  must  be  represented  in  a  somewhat  different  form. 

Differentiating  I  for  this  purpose  2s  times  with  respect  to  t  and  consid- 
8 

ering  (12) ,  we  obtain 

<*?•  “  R 


where 


As  is  easy  to  check  by  differentiation,  the  last  equation  is 
satisfied  by  the  function 

(2,  —  l)|  i  ~R~  ^ “  *)’*  * & coJ  $  - - -  p 


where  R ^  =  Jr?  +  ? . 


(a)  (a) 

For  »  C1  *  • • • »  by  sequentially  differentiating  the  last 

expression  and  assuming  £  ■  0  we  obtain 


and  in  general 


(«i _  1  (dl  f,  \  _ (—  O'  f 

c'l)— »[!?),=  -  /«.  1 


Using  the  formulas  given  by  Watson  £10]  for  Integrals  of  this  type 
(see  also  the  paper  by  Fock  [11]),  we  obtain 


C‘V>==  2  {w) 


(—  l)'c' 1 H  ,_i(feo>) 


nr 


(-¥) 


-  li  - 


The  variables.  A  and  D  used  in  the  text  are  equal  respectively  to 
c0^and  By  replacing  the  Hankel  functions  by  their  asymptotic 

representations,  from  the  last  formula  we  obtain  the  values  of  (15)  for 
them. 
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